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The visible spectrum of is studied using high-sensitivity action spectroscopy in a cryogenic radiofrequency 
multipole trap. Advances are made to measure the weak ro-vibrational transitions from the lowest rotational 
states of up to high excitation energies providing visible line intensities and, after normalisation to an 
infrared calibration line, the corresponding Einstein B coefficients. Ab initio predictions for the Einstein B 
coefficients are obtained from a highly precise dipole moment surface of and found to be in excellent 
agreement, even in the region where states have been classified as chaotic. 


The triatomic molecular ion prevails in cold hydro¬ 
gen plasmas and has therefore been the subject of study 
for many decades in astronomy and for even longer in 
molecular physics. In space, H;/ is the cornerstone of 
the chemical network leading to larger moleculeJi^ and it 
is now being widely monitored in both the dens^ and 
diffus^ interstellar medium. In planetary atmospheres, 
it acts as a tracer of H 2 , as an important coolant, and as 
a tracer of ionospheric activity.l^ H;/ is thought to be the 
key of the thermal stability of hot Jupiter exoplanets.^ 
For modeling the role of H;/ in these environments and 
for interpreting related observations, precise and reliable 
data on the H;/ spectrum are required in an increas¬ 
ingly wide wavelength range. Such work relies heavily 
on comprehensive, theoretical spectral line calculati ons 
including transition energies and also line intensities.*^ 
While theoretical line intensities in comparison to labo¬ 
ratory observations served as a guide in the assignment 
of complex discharge spectra in the infrared (IR) spectral 
region,*^ experimental information on the intensities for 
independently identified lines is scarce. 

By combining laser spectroscopy in a cryogenic trap 
with high-precision calculations of the ground-state po¬ 
tential energy surface of the ion, we recently identified 
several lines of in the near-visible and visible (VIS) 
spectral range and showed that our theoretical transition 
energies agreed with the measured line positions on the 
sub-tenth cm“^ level.l^Here we prepare H;/ ions in their 
lowest ortho and para rotational levels only and measure 
the intensities of eighteen weak near-VIS and VIS spec¬ 
tral lines, which are directly referred to the intensity of 
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an IR calibration transition whose value can be safely 
predicted. Theoretically, we developed a dipole moment 
surface (DMS) for H;/ that also covers the VIS part of 
its vibrational spectrum, which we used to model the line 
intensity measurements. 

Highly excited molecular motion forms a challenge 
in the ab in itio t reatment of molecules and molecu¬ 
lar dynamics! ^*^ * ^^ * The Born-Oppenheimer approxima¬ 
tion breaks dowrfi^ and the molecular motions couple 
strongly to each other, leading to ill-defined approximate 
quantum numbers and vibrational angular momenta.l^ 
This behavior is associated with the onset of chaos in 
the analogous classical problem. Non-adiabatic effects 
are often only included as ad hoc correction factors on 
the calculations.!^ Despite these difficulties, the recent 
advanced theoretical treatmentP gave spectroscopic accu¬ 
racy in the predictions of ro-vibrational energy levels 
up to 2 eV internal excitation. 

The investigation of high levels of excitation also con¬ 
stitutes a great experimental challenge as the spectra 
often beco me complex and unassignable. Carrington 
et al^^^ triggered the extensive theoretical investiga¬ 
tion of with their measured, rich spectrum of states 
near the dissociation limit, which remains unassigned. 
The H;/ internal motion already becomes chaotic at 
about one-third of the dissociation energy,!^! where the 
nuclei sample linear configurations. The first measure¬ 
ments just above this barrier to linearity were performed 
by Gottfried and coworkers using advanced absorption 
spectroscopy.*^^^*^ This method reaches its limit at 13 700 
cm“ ^, where the absorption Einstein coefficient B of the 
transitions reaches a few 10^® cm® J“^s“^. No Einstein 
coefficients could be derived due to unknown initial pop¬ 
ulation of the H;/ levels. The identification of the many 
lines was also hindered by the large number of initial 
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TABLE I. Observed near-VIS and VIS transitions of with their assignments, measured transition wave number Vm, and 
their intensity relative to the IR reference line. The transitions start from the vibrational ground state OO'^. 


AJ(J",G") 


Vruicra ^) Intensity ratio 

This experiment 

B (10^® cm® J-%-®) 
This exp. p| This theory 

NMTpl 

R(1,0) 

22® (2,3) 

10 752.161(5) 

0.047(10) 

73(16) 

75.9(78) 

60.5 

^^(1,1) 

05^ (0,1) 

10 798.777(5) 

0.0172(19) 

26.5(30) 

34.6(12) 

32.9 

g(i,0) 

05^ (1,0) 

10831.681(5) 

0.073(10) 

112(16) 

98.2(37) 

94.0 

Q(i,i) 

23^ (1,1) 

12 373.526(10) 

0.002 77(64) 

4.3(10) 

4.87(83) 

4.04 

R(l,l) 

23^(2,1) 

12 381.137(10) 

0.002 68(66) 

4.1(10) 

4.67(52) 

4.10 

^^(1,1) 

06® (0,2) 

12 413.247(10) 

0.002 98(75) 

4.6(12) 

3.85(28) 

3.73 

R(l,l) 

23^(2,1) 

12 588.951(5) 

0.000 72(25) 

1.10(38) 

1.05(18) 

0.86 

R(l,l) 

23® (2, 5) 

12 620.223(5) 

0.004 07(80) 

6.3(12) 

5.09(31) 

4.86 

R(l,l) 

06® (2,4) 

12 678.683(10) 

0.005 5(11) 

8.6(17) 

8.74(49) 

8.01 

R(l,l) 

32® (2,4) 

13 332.884(10) 

0.002 61(84) 

4.0(13) 

2.86(40) 

2.05 

^’(1,1) 

07^ (0,1) 

13638.251(5) 

0.002 54(95) 

3.9(15) 

4.47(44) 

4.14 

^^(1,1) 

08® (0,2) 

15 058.680(5) 

0.00100(22) 

1.53(33) 

1.76(27) 

1.62 

Q(l,l) 

16® (1,2) 

15 130.480(5) 

0.000 47(10) 

0.72(16) 

0.83(22) 

0.71 

R{1,0) 

16® (2, 3) 

15 450.112(5) 

0.000 489(66) 

0.75(10) 

0.92(31) 

0.78 

R{1,0) 

16^ (2, 3) 

15 643.052(5) 

0.000 718(98) 

1.11(15) 

1.18(30) 

1.01 

Q(l,l) 

25^ (1,1) 

15 716.813(5) 

0.00104(33) 

1.60(51) 

1.56(24) 

1.40 

R(1,0) 

34® (2,3) 

16 506.139(5) 

0.000 83(32) 

1.28(50) 

1.23(23) 


Q(l,l) 

J' = lfl 

16 660.240(5) 

0.000 24(12) 

0.38(19) 

0.62(18) 



^ Using B = 1.54(2) X 10^^ cm® J ® s ^ (this theory) for the reference line 03^ <— 00® i?(l, 1)* at 7144.235 cm ^ 
® Reference [7| 

Further assignment unavailable 


levels involved. 

The development of action spectroscopy in a cryogenic 
radiofrequency (RF) ion trap provided measurements in 
the near-IR range^^ and showed great promise for attain¬ 
ing much higher sensitivity than absorption spectroscopy. 
We have recently extended the range of this method up 
to almost 17000 cm“^.l^ The measurement of line inten¬ 
sities in the VIS region has, howev er, bee n out of reach 
until now. As described previousljEEsHni ions were 
stored in a multipole RF ion trap and buffer-gas cooled 
with He (2 x 10^^ cm“^) to ~ 85 K. About 85% of the 
ions thus populated the two lowest rotational (J = 1) 
levels of the vibrational ground state, having ortho and 
para symmetry. The near-VIS to VIS transitions were 
measured by chemical probing with Ar (1 x 10^^ cm“^). 
The reaction -|- Ar —>■ ArH+ -|- H 2 is endothermic 
but activated by ro-vibrational excitation of above 
4500 cm“^. The count of extracted ArH+ reaction prod¬ 
ucts per storage cycle is the spectroscopic signal and is 
only formed by ions excited above the endothermic- 
ity barrier. There is a background from non-laser in¬ 
duced ArH'*', originating from the initially hot H;^ ions. 
The decay rate of this background ArH’*' is determined 
by the back reaction of ArH+ with H 2 (2 x 10^° cm“^), 
while the forward reaction is stopped within a few ms 
by the buffer-gas cooling of H;^. A pre-cooling time of 
200 ms was chosen to achieve a sufficiently low ArH+ 
background level of less than 0.1 counts on average per 
storage cycle. Afterwards, the ions were irradiated 
for 100 ms yielding, when on resonance, internal ion 


energies much higher than the required reaction barrier. 
The reaction time of ^0.5 ms for ArH+ formation, esti¬ 
mated from the Langevin cross section, is shorter than 
the effective radiative decay time of at least several ms.^ 
Collisional de-excitation of vibrations can be neglected, 
even for He which has the highest partial pressure. After 
irradiation, the ion cloud was released into a mass spec¬ 
trometer (MS) that could be set to either ArH+ and H;)" 
detection from one 300-ms storage cycle to the next. The 
ions were counted on a Daly stamp with a scintillator and 
a photomultiplier detection systenP^HH that is sensitive 
enough to single-ion count the low ArH'*' and the high 
H;^ (normalization) ion numbers. 

To determine the intensity of the spectral lines, we 
performed scans of the laser wavelength in small steps 
around the central wavelength pre-located in wider 
search scans.l^ The strength of the ArH+ counting sig¬ 
nal exceeding the background level was determined from 
a Gaussian fit to the scanned line. The observed lines 
were Doppler-broadened corresponding to a kinetic tem¬ 
perature near 85 K. The experimental challenge for reli¬ 
able intensity measurements is not only to increase the 
spectroscopic sensitivity to a level where we can measure 
these very weak, line^but also to increase the dynamical 
counting range for both the ArH+ signal and the H;^ nor¬ 
malization signal, whilst retaining linearity and constant 
(within the Poisson noise) experimental parameters. To 
increase the dynamic range to 10"*^, we applied a volt¬ 
age ramp on the trap exit electrode to gradually release 
the ions, leading to a stretched ion pulse train of 0.5-1 
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ms and yielding a factor of 25 improvement in dynamic 
range. Group recordings of 10 storage cycles each were 
alternately performed with the MS set on (normaliza¬ 
tion) or on ArH+ with (signal) and without (background) 
laser. Two spectroscopy lasers were applied along the 
trap axis: a single-mode Ti:sapphire or a dye laser for 
measuring a near-VIS or VIS line, and a single-mode IR 
laser diode for measuring an IR reference line, 03^ ^ 00° 
i?(l,l)* at 7144.235 cm“^.l^ Both lasers filled identical 
apertures, ensuring their spatial overlap inside the trap¬ 
ping region. A measurement on a near-VIS or VIS line 
typically lasted 15 h, and immediately before and after 
this, 2-h measurements were performed on the IR refer¬ 
ence line. In a campaign extending over several months, 
intensity measurements relative to the IR reference line 
were performed on 18 near-VIS or VIS lines of These 

intensity measurements were performed separately from 
the runs aiming at the line frequency determination de¬ 
scribed earlier.El 

The intensity ratios, linking the near-VIS to VIS lines 
to a single IR calibration line far below the barrier to lin¬ 
earity, are listed together with the measured transi tion 
wave numbers Vm in Table |l] In the line assignments ,121221 
J" and J' refer to the lower and upper level, respectively. 
Also listed are values of the absolute rate coefhcients ob¬ 
tained for the Einstein coefficient of the IR calibration 
line from the theoretical method of the present work; pre¬ 
vious theoretical calculation^^ of this calibration value 
yielded essentially the same result {B = 1.55 x 10^^ 
cm° s“^). 

In our calculations on the two-electron system the 
spatial component of the electronic wave function $M(r), 
is expanded in terms of a set of M explicitly correlated 
Gaussians (EGGs), {gk}k=i,...M, as: 

M 

^’M(r) = ^ CkVgkir), (I) 

k=l 

where V is the permutational-symmetry projection op¬ 
erator which symmetrizes each EGG with respect to the 
permutation of the electron labels. The EGGs are the 
following 6 -dimensional functions: 

5 fc(r) = exp [-(r - Sfc)'Afe(r - Sfc)] , (2) 

where r and Sk are 6 -dimensional vectors of the elec¬ 
tronic Cartesian coordinates and of the coordinates of 
the Gaussian shifts, respectively, and the prime denotes 
vector transposition. Aj, is a 6 x 6 symmetric matrix of 
the Gaussian exponential parameters defined as: 

Afc = Afc 0 I 3 , (3) 

with I 3 being the 3x3 identity matrix and 0 denoting 
the Kronecker product and Ak is a 2 x 2 symmetric ma¬ 
trix. More details on how the spatial wave function is 
constrained to satisfy the Pauli exclusion princi ple, and 
to be square integrable are given elsewhere 



Transition wave number (cm-') 


FIG. 1. Line intensities for near-VIS and VIS transitions of 
H 3 + as a function of the transition wave number (horizontally 
shifted where helpful for clarity). Experimental results (dots 
and error bars) derived from the measured line intensity ratios 
(right-hand scale) and theoretical results from this work (open 
squares) and from NMT® (full squares, where available). 

Here calculations of the dipole moment surface (DMS) 
are presented. The dipole moment is determined as an 
expectation value of the dipole operator with the varia- 
tionally optimized electronic wave function 0 contai ning 
900 EGGs for each of roughly 42 000 geometry points.l^I^ 
The dipole moment vector associated with a certain 
molecular configuration defined by the Cartesian position 
coordinates (Ri, R 2 , and R 3 ) is: 

M / 2 \ 3 

ti = '^CkCi{Vgk{r)\ \Vgiir)) +'^Ra 

kl \ i=l / a=l 

M / 2 \ 3 

= ^CfcQ(g,(r)| ^(r)) + ^R„,(4) 

kl \ i=l / a=l 

where are the coordinates of electron i, R„ are the 
coordinates of nucleus a both expressed in a Cartesian 
coordinate system whose origin is the center of the posi¬ 
tive charge of the H 3 system. In eq. 0 we use the fact 
that the electronic dipole-moment operator, — 
is symmetric with respect to the permutation of the elec¬ 
trons and V = W ■ The spatial wave function is assumed 
to be normalized. The dipole moment matrix element in¬ 
volving two EGGs, 5fc(r) and giiy) is given asP^ 

{gk{r)\ri\gi{r)) = S'fcij [A^ -f A;] ^ [A^s^ -b A/S;] | , 

(5) 

where {v}i denotes the f-th component of vector v. Ski = 
{gk\gi) is the overlap integral. 

Since the fits to produce a DMS were not entirely stable 
we only used dipole components from 4257 points lying 
up to 30 000 cm“^ above the minimum generated by the 
EGG calculations in the form and /ij,. Previous 
DMSs were based on a grid of only 69 points; following 
these studiesthe dipole moments were transferred 
to effective charges dx and dy which were expanded into 
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a power series in terms of symmetry coordinates. Since 
our DMS showed sensitivity to the htting procedure, we 
opted for the simplest form and used only 7 constants 
and the function of Ref. nzi which reproduce the com¬ 
puted dipoles for points lying below 18 000 cm“^ with 
a standard deviation of 0.001 a.u. Appropriate func¬ 
tional forms for fits of a global DMS will be the subject 
of further work. Transition intensities were computed us¬ 
ing nuclear motion program DVRSthe non-adiabatic 
model of Polyansky and TennysorP^ with other param¬ 
eters chosen as i n our previous studies of Hg VIS tran¬ 
sition frequencies .12*23] major uncertainty in the 

calculations is given by the fit of the DMS, we used the 
extremes of plausible fits to estimate uncertainties in our 
calculated intensities; these uncertainties (see Table 1) 
are systematic rather than statistical errors and roughly 
represent a l-u conhdence level. 

The high agreement between experiment and the 
present calculations is shown in Table |T] and Fig. and 
covers the entire range of B coefficients whose magnitude 
varies by a factor of 300. For comparison, the feure also 
shows the results of Neale, Miller and TennysorP (NMT), 
the only previous intensity calculation extending into the 
VIS. The NMT results also agree well with the observa¬ 
tions which is somewhat surprising given that the study 
was not designed to be accurate above the barrier to lin¬ 
earity and that NMT’s frequencies are signihcantly less 
accurate than the ECG ones.P 

spectra provide a unique probe of the Universe but 
their interpretation relies heavily on the reliability of ab 
initio intensity calculations. The present studies suggest 
that, with care, the intensities can be computed over a 
wide wavelength range. Following the confirmation of 
the accuracy of our ab initio model for both tran¬ 
sition frequencies and, with the present study, intensi¬ 
ties, work on a new, experimentally tested, comprehen¬ 
sive line list for astronomical and other applications 
is now in progress. This will make highly excited lev¬ 
els of H;^ broadly available as a tool of precision laser 
spectroscopy, enabling ultra-sensitive spectroscopic de¬ 
tection of by multistep photodissociation.l^ It will 
also pave the way to the reliable modeling of broadband 
spectral emission from in astrophysical environments, 
as considered regarding the formatioiP^ and evolutiorl^ 
of metal-poor stars, and the stability of the atmospheres 
of hot exoplanets.l^ Furthermore, the use of accurate in¬ 
tensities for interpreting the Carrington-Kennedy near¬ 
dissociation will allow the division of the spectra in 
hierarchies, greatly simplifying its analysis.l^ 
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